Gain-assisted slow to superluminal group velocity manipulation in nano-waveguides 



o 
o 

< 

(N 



O 

• i-H 

-f— > 

Or 

q . 

o ■ 

>v 
43 ■ 

Oh: 



> ■ 
■ 

o\ : 
o . 

00 ■ 

o ■ 
\o : 
o . 

o : 
43 : 

Ph. 



X 



Alexander A. Govyadinov and Viktor A. PodolskijQ 

Physics Department, Oregon State University, 301 Weniger Hall, Corvallis, OR 97331, USA 

We study the energy propagation in subwavelength waveguides and demonstrate that the mecha- 
nism of material gain, previously suggested for loss compensation, is also a powerful tool to manipu- 
late dispersion and propagation characteristics of electromagnetic pulses at the nanoscale. We show 
theoretically that the group velocity in lossy nano-waveguides can be controlled from slow to super- 
luminal values by the material gain and waveguide geometry and develop an analytical description of 
the relevant physics. We utilize the developed formalism to show that gain-assisted dispersion man- 
agement can be used to control the transition between "photonic-funnel" and "photonic-compressor" 
regimes in tapered nano-waveguides. The phenomenon of strong modulation of group velocity in 
subwavelength structures can be realized in waveguides with different geometries, and is present for 
both volume and surface-modes. 



Conventional optical fibers support propagating modes 
only when waveguide radius is sufficiently large [lj. In 
contrast to this behavior, plasmonic systems, anisotropy- 
based waveguides, nanoparticle chains, and optical coax- 
ial cables 3, S S 13 support energy propagation even 
when the typical waveguide cross-section is much smaller 
than the wavelength. Unfortunately, since the major- 
ity of these nano-waveguides rely on plasmonic materials 
to confine the radiation beyond the diffraction limit, the 
propagation of nano-constrained radiation is often lim- 
ited by material losses. While the emerging field of ac- 
tive plasmonics |(| promises to overcome absorption lim- 
itations in nano-waveguides, full compensation of losses 
appears to be experimentally challenging |2j. 

In this Letter we focus on gain-assisted phenomena 
beyond absorption compensation and study the perspec- 
tives of controlling the dispersive properties of active 
nanoscale waveguides. We show that even relatively weak 
material gain, which is unable to compensate losses, is 
capable of producing large variations of the group ve- 
locity, bringing such exotic phenomena as slow (0 < 
v g < c) and ultra-fast (v g < 0) lightjl H E3 to the 
nanoscale domain. However, in contrast to diffraction- 
limited systems, where the group velocity is controlled 
solely by material dispersion, the energy propagation in 
nano-waveguides is also strongly affected by the waveg- 
uide geometry. We demonstrate that interplay between 
geometry- and material-controlled modal dispersion in 
tapered waveguides leads to the transition between pho- 
tonic compressor regime, where the reduction of phase 
velocity is accompanied by the simultaneous reduction 
of group velocity [111 Il2| and photonic funnel [l3| regime 
where the product of phase and group velocities re- 
mains constant . The developed formalism is illustrated 
on the examples of two fundamentally different nano- 
waveguides: a surface-mode-based plasmonic nanorod, 
and a volume-mode fiber with anisotropic core. Appli- 
cations include nanosized tunable delay lines, all-optical 
buffers and data synchronizers. 

The geometry of a typical waveguide structure is 
schematically shown in Fig. For simplicity, we as- 



sume that the system is homogeneous along the propa- 
gation direction (z axis) the transverse structure of 
the waveguide is represented as a set of N layers with 
different permittivities. Any wavepacket, propagating in 
such a structure can be described as a linear combination 
of waveguide modes, with each such mode having its own 
dispersion relation between propagation constant k z and 
frequency ui. 

To find this relation, we first represent the field in 
i-th layer by a linear combination of "incoming" and 
"outgoing" (in radial direction) waves with the same z- 
dependence, but different structure across the layer (de- 
scribed by parameters Ki)0- We then relate the fields in 
each layer to the fields in the neighboring layers through 
the boundary conditions, leading to the waveguide ma- 
trix M that relates the (complex) amplitudes of waves 
in all layers. Since a waveguide mode corresponds to a 
nontrivial solution of the system of algebraic equations 
described by M, the modal dispersion equation can be 







cladding -r 


"V x 






core k * 











0.55 A,|jm 0.075 
0.025 





FIG. I: (color online) (a) The schematic geometry of a mul- 
tilayered waveguide; two-layer cylindrical system with axial 
symmetry is shown, (b) Dielectric permittivity of Ag; real 
and imaginary parts are shown in solid and dashed respec- 
tively, (c) Real (top) and imaginary (bottom) parts of the di- 
electric permittivity of Rhodamine-6G model (see text). Solid 
blue, dashed green, dash-dotted brown and dash-dot-dotted 
red curves correspond to gain values 0%, 33%, 66% and 100% 
respectively. 
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found using: 



i(k z ,UJ,K) 



detl 



0. 



(1) 



where k is a vector constructed of Ki, which in turn are 
found from the wave- and layer-specific dispersion equa- 
tions for individual waves: 



A(fc z ,w,K t ) = 0, 



(2) 



Eas. (|ll2|) provide a complete description of electromag- 
netism in multi-layered waveguides. In particular, the 
group velocity is given by the following expression: 



&B_ ( dDj \ 1 dDj _ 3D 
dni \ dKi J dk z dk z 

dP ( dDj \ ~ 1 dDj dB 
dhii \ dni J duj duj 



(3) 



The number of terms in the summations in Eq. @ (the 
length of k) depends on the waveguide symmetry and 
individual layer properties. Thus, when the system has 
axial symmetry (see Fig. ^i), each layer with isotropic 
permittivity adds one k per "axial number" m, each layer 
with uniaxial anisotropy adds two such terms, etc. 

Eq. © can be used for an arbitrary waveguide ge- 
ometry including planar, square, and circular systems, 
and can be applied to plasmonic, coaxial, and volume 
waveguide modes. To illustrate the developed formal- 
ism we consider two axially-symnietric cylindrical nano- 
waveguides: (I) a plasmonic nanorod in the dielectric 
material, and (II) a waveguide with anisotropic core, also 
known as photonic funnel |l3j . 

In the case of the plasmonic waveguide the z 
component of the m-th cylindrical wave inside and 
outside the metal nanorod can be represented by 
I m ( Kl r)e im 't' +lk - z -^ t and K m {K 2 r)e lm< t> +lk > z - lut respec- 
tively, with I m and K m being modified Bessel functions, 
and k\ 2 = kl — e\^ 2 1 ' <? . The boundary-matching tech- 
nique described above yields the following set of disper- 
sion equations for the first (m = 0) TM surface plasmon 
polariton (SPP) mode 0: 



nSPP 



d spp 
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OJ n 2 



(4) 



Direct substitution of Eqs.(@J into Eq.© yields: 
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with v p being phase velocity of the mode: 
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FIG. 2: (color online) Group and phase velocities of an SPP 
on metal nanorod in Rh6G methanol solution (see text) as 
functions of gain, frequency (a,c) and radius (b,d). The radial 
dependence is given at A = 534nm, the wavelength one - at 
R = 35nm. 



When the nanorod radius is much smaller than the 
free-space wavelength (R <C Ao), K1.2 is proportional to 
1/R. In this regime the phase velocity of SPP is propor- 
tional to the radius, the terms involving D SPP in Eq. (0 
are proportional to R 3 , while the remaining term is pro- 
portional to R. 

The dramatic effect of material dispersion on the re- 
lationship between the phase and group velocities is 
now clearly seen. In the absence of material dispersion, 



Jg vp 



const, and the decrease of phase velocity with 
the radius is accompanied by the increase of group ve- 
locity. This effect, introduced for anisotropic waveguides 
(see below') [1 3| is somewhat similar to the motion of un- 
compressible fluid through a funnel: the decrease of the 
cross-section leads to the increase of the local speed . 

Any non- vanishing dispersion, however, dramatically 
changes this behavior and leads to the regime when both 
v p and v g simultaneously vanish. This regime, that origi- 
nates from the material dispersion of plasmonic compos- 
ites, has been shown to yield slow light in nanoscale plas- 
monic wave guid es | ll| and in adiabatic plasmonic energy 
compressors [l2j . 

The crucial point of our work is that the dispersion 
of the dielectric core can be used to control the group 
velocity of SPP by adjusting the term — i^^J 

in Eq. This way the nanoplasmonic system switches 
between "photonic funnel" and "photonic compressor" 
regimes. Note that the group velocity can be indepen- 
dently controlled by either material dispersion or waveg- 
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uide radius. It is therefore possible to build plasmonic 
systems with either "fast" or "slow" modes (see Fig|2J, 
or implement an adjustable gain mechanism to tune in 
between these two regimes. It is also possible to con- 
struct a tapered plasmonic fiber, similar to that described 
in|l2|. in which the packets would travel with superlumi- 
nal speed at larger radii, and compress toward the small- 
radius apex. 

To further demonstrate this behavior and show the 
control over the group velocity with the material gain 
we model the dispersive properties of a silver nanorod 
(described by Drude model, Fig^>): 



ex(uj) 



As 



Ab 



(7) 



submerged into the 10% solution (0.1-M) of Rhodaminc- 
6G in Methanol (Rh6G), which in the optical frequency 
range can be approximated by: 



62 M 
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Rh6G , 
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In our calculations we use e^? 



5, UJ p - 
,Rh6G 



46.26/zm- 1 , 

7 P = 0.11/im _i for silver [18| and e™ ou = 1.81, 71 = 
0.4/xm -1 , 07( ) 1 = 12.82iitm~ 1 , 72 = 0.2^m _1 , W( ) a = 
11.74/im -1 for Rhodamine |l9j. We model the material 
gain by fixing A\ — 0.001 and adjusting A 2 to gradually 
change the corresponding resonance strength. On the 
microscopic level, this process corresponds to a gradual 
increase of the population of the excited level of Rh6G 
with respect to the ground-state population, achievable, 
for example, by an external pump. The dependence of 
the dielectric permittivity of Rh6G on the material gain 
measured in percents of the excited state population is 
illustrated in Fig. In the remaining of this work, 
we assume the regime with less than 50% of population 
in the excited state, corresponding to A 2 > (A 2 < 
represents an inverted system). 

The group and phase velocities of SPP at the Ag-Rh6G 
interface are shown in Fig. [21 As mentioned above, the 
gain level used in our work is relatively weak. Indeed, 
the phase velocity is almost unaffected by the material 
gain (see Fig.|2t,d). Furthermore, the imaginary part of 
62, and of propagation constant k" remain positive, in- 
dicating that the gain is insufficient to compensate SPP 
losses. Both superluminal (v g < 0) and slow (0 < v g -C c) 
light regimes can be identified in active nanoplasmonic 
structures. We note that while the SPP mode remains 
lossy, the total absorption in the fast- and slow-light re- 
gions is relatively small. The corresponding parameter 
\K/k z \ < 0.05. 

When the gain is close to 25%, the group velocity un- 
dergoes the transition from superluminal to slow light at 
a constant pumping level with varying nanorod radius. 
Similarly, at R — 35nm, the superluminal to slow light 




FIG. 3: (color online) Group and phase velocities of the TM01 
mode in anisotropy-based cylindrical subdiffraction waveg- 
uide (see text) as functions of gain, frequency (a,c) and radius 
(b,d). Radial dependence is given at A = 534nm, the wave- 
length one - at R = 35nm. Note that the phase velocity is 
negative. 



transition happens as a function of pumping, providing 
almost unlimited control over the group velocity of SPPs 
at the nanoscale. 

To demonstrate universality of nanoscale group ve- 
locity modulation with respect to waveguide modes, 
we apply the developed formalism to another type of 
subdiffraction waveguide - the photonic funnel. These 
anisotropy-based systems support volume modes with ei- 
ther positive or negative refractive index |2(lj and can be 
easily integrated with diffraction-limited fibers^j|. I n 
the simplest case of a TMqi mode in cylindrical geome- 
try with uniaxial dielectric core with optical axis parallel 
to the direction of mode propagation and perfectly con- 
ducting metallic walls, Eqs. (1112ft become: 



— Jo(kiR) 

1, 1 2 k 2 k 2 



(9) 



where e xy and e 



are the components of effective dielectric 
permittivity along and perpendicular to z-direction (see 
Fig.QJi and ref. [3J) and Jo is a Bessel function of the first 
kind. In the nanoscale limit (R <C Ao) Eq. becomes: 



VnV v 
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1 de z 
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(10) 



Similar to what has been shown earlier for SPP mode, 
this expression represents the interplay between the "fun- 
nel limit" VgVp = c 2 /e xy = const and the "compressor 
limit" VgVp (x l/«i oc R 2 . 
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Fig. |31 shows the behavior of the TMqi mode in the 
waveguide with metamaterial core, consisting of alter- 
nating (in z direction) layers of Ag and Rh6G. In our 
computations we use Eas. (l7l8l) to calculate the dielectric 
constants of Ag and Rh6G, accompanied by the effective- 
medium technique to obtain the effective permittivities 
of the nanolayer material (see Ref. 3] for details). 

Note that the phase velocity (as well as the effective 
refractive index) of the metamaterial photonic funnel is 
negative 0| . Apart from the sign of v p , the propaga- 
tion of volume mode in metamaterial structure (Fig. I3J 
is similar to that of an SPPs (Fig. 0) . Once again, one 
can effectively control the group velocity between ultra- 
fast (v g < 0) and slow (0 < v g <C c) values in lossy 
(\k"/k z \ < 0.08) nanoscale system by tuning the mate- 
rial dispersion of the core or waveguide radius. 

To conclude, in this Letter we demonstrated the pos- 
sibility of versatile (from slow to superluminal) group 
velocity modulation in sub-diffraction waveguides. The 
developed formalism can be directly utilized in different 
waveguide geometries (planar, square, etc.), different set 
of materials (polar dielectrics, semiconductors, quantum 
wells [2lj. quantum dots H2), and scaled for UV, opti- 
cal, IR, THz and microwave spectral regions. Finally, we 
note that the strong modulation of group velocity can be 
also achieved in fully-active (amplifying) media, where 
it will be accompanied by loss-less or amplifying mode 
propagation^. 
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